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INFLUENCE OF MAGNETIC SHEAR ON THE COLLISIONAL CURRENT DRIVEN 
ION CYCLOTRON INSTABILITY 


I. INTRODUCTION 

Intense field aligned currents observed In the auroral G and F regions 
(Park and Cloutier, 1971, Cahill et al., 1974, Kelley et al., 1975) were 
conjectured to be responsible for the generation of electrostatic ion 
cyclotron waves. Since the ionosphere is collisional in the E and F regions 
the field aligned currents give rise to negative energy waves which grow due 
to the dissipation associated with electron-neutral collisions (Chaturvedi, 
1976). These are electrostatic waves and lead to density irregularities in 
the medium. In fact, radar measurements in the E region detect 

irregularities collocated with auroral electrojets (Greenwald et al., 1975), 
and a.c. field measurements (Kelley et al., 1975) have been identified as 
being due to unstable ion cyclotron waves (Drummond and Rosenbluth, 1962 - 
collisionless domain; Kindel and Kennel, 1971 - weakly collisional 

domain). Sounding rockets launched from Syowa station, detected 
transversely propagating electrostatic plasma Instabilities of scale sizes 
around 200 km in association with strong field aligned currents (Ogawa et 
al., 1981). Chaturvedi (1976), using a local fluid analysis, showed that, 
in a partially ionized plasma, field aligned currents can support almost 
transversely propagating ion cyclotron waves owing to resistivity 
experienced by electron parallel motion (electron-neutral collisions), while 
the ion-neutral collisions were found to have a mild stabilizing effect. 
The field aligned currents observed are usually '■ 1-10 pA/m 2 (Kelley et al., 
1975), although recently Burke et al. (1983) reported observations of large 
2 

currents ** 135 pA/ra . These field aligned currents can generate magnetic 
shear which in turn affect the mode structure. This is especially so in the 
case of collisionless current driven ion cyclotron waves (Gangull and 
Bakshi, 1982). 

In this brief report we investigate the effects of a self-consistent 
magnetic shear on the collisional current driven ion cyclotron instability 
(CIC). We find that the magnetic shear has no noticable stabilizing effects 
on the current driven ion cyclotron instability in the domains of interest 
pertaining to auroral conditions (characterized by V, ' 0.5 - 5 km/s, 

o a 4-1 

corresponding to currents of order 10 - 100 pA/tn , and v @ -*• 10 s , where 
Manuscript approved March 20,1984. 
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is the electron drift velocity along the magnetic field and is the 

electron-neutral collision frequency). This is because the self-consistent 

magnetic shear produced by the currents observed in the auroral region is 

very small. This result is in contrast to the major effects found by 

Ganguli and Bakshi (1982) in the collisionless domain. Their kinetic 

treatment shows that, due to the nonlocal boundary conditions, magnetic 

shear affects the collisionless current driven ion cyclotron instability in 

two ways: One, due to an explicit shear dependent term of 0(p g /L g ) which 

vanishes for P g /L g + 0* Two, due to a shear independent term containing the 

derivatives of the potential function defined in equation 18. In the fluid 

example, the second term happens to be zero (since the derivative of Q, 

defined in Eq. 18, vanishes at Xq which is defined in Eq. 19) and hence the 

lack of the significant effects in the p /L ♦ 0 limit. Furthermore, 

s s 

for P g /L g « 1 magnetic shear localizes the wave packet around k^/k « 1, 
indicating that the mode is almost perpendicularly propagating, and stronger 
magnetic shear stabilizes the instability. The paper is organized as 
follows. In the next section we give the theory and derive the nonlocal 

mode structure equation. In the third section we present the results, and in 
the last section we apply the results to the auroral ionosphere and discuss 
future work. 

II. THEORY 

We consider a partially ionized plasma such as the one encountered in 
the auroral E region. The geometry used in this analysis is as follows. 
The magnetic field is aligned with the z-direction along which flow the 
currents that drive the CIC. Because of the self-consistent shear produced 
by the field aligned currents, the magnetic field acquires a component along 
the y-direction as x is varied, 

B(x) » B {z_ + x/L g j). (1) 

This implies that the parallel wave vector becomes a function of x. 


k =* 


k z + V x/ V’ 


( 2 ) 







where L g is the shear length. The self-consistent shear is calculated from 
the Maxwell's equation 


7 x B - — J - — (.neV.J (3) 

— c — c —a 

where is the drift velocity parallel to the z-direction, and is given by 
L s “ a ’ lp s^ c s /V d^ <4) 


a s U 8 /cJ 2 <M/mHa. 2 e /n 2 J 


(5) 


where <»> and fl are the electron plasma and cyclotron frequencies, 

pe e 

respectively, c is the speed of light, and M and m are the ion (NO ) and 
electron masses, respectively. 

The basic equations describing the problem are as follows. The 
electron and ion continuity equations are given by 


ir + ■ 0 (6 > 

where a indicates the species. The equation of motion for the ions is. 


‘ 3F h m ' ’"i + en l- 7 * + > - ltov i h- 


(7) 


while the equation of motion for electrons is, 

V x B 

0 “ - Vp e - en(-V$ + ~~5 - —*J - mnv^ , (8) 

where \> @ and are electron-neutral and ion-neutral collision frequencies, 
respectively, p a » n a ^ a * s c ^ e pressure of species a, and d/dt * 

3/3t + V^V. In this paper we ignore the electron inertia, and assume that 

« 1* and v^/0^ « 1. We also ignore the viscosity and thermal 
conductivity terms because these terms only introduce a numerical constant 
that multiplies the growth rate and do not introduce any additional k 2 /ky 
terms (see, for example, Chaturvedl and Kaw, 1975). We also ignore the 
elctromagnetlc effects and consider only to electrostatic perturbations. We 
derive the mode structure equation assuming that the perturbed quantities 
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vary as f ' f(x) exp |.-i(o»t - k y y - k^zJJ . We note that we do not Fourier 
analyze in the x direction since the magnetic field is assumed to vary in 
that direction. 

From Eq. (8), we obtain the perpendicular and parallel (to B) 
components of the electron velocity 

- (_i_5L . £ 7 <b 1 _®| 

1 mn m ft J 


i v P 

. _ e . 

v 1 mn m z* J 


From Eq. (10) we define 


d “ ez v '-B zT 


From Eq. (7) we obtain for the ions 


^■[i + (i.»V 2 /a*J U ‘ V + “i 


where we have assumed T g **- , and <f> - $° + <p (<j> being the perturbed 
potential). Substituting Eqs. (9) thru (13) in Eq. (6) we obtain, after 
subtracting the electron continuity equation from the ion continuity 
equation, the nonlocal mode structure equation 


2 2 2 r m 

3x k z c 8 - Hm/Mjv e (u-k z V d J 


local limit of Eq. (14) (L •* - and 3/3x « k ) yields Q(x) - 0, where Q(x) 







are the second two terms in the above equation. This gives the results 
obtained previously by Chaturvedi (1976). The linear local growth rate is 


k k • V, 

■ 2-M - =-=1) v 

.2 1 (D 1 @ 


rn 2 . n 2 + k 2 C 2 . 
r i y s 


Gq. (15) leads to the threshold condition 


Note that the nonlocal treatment allows for the formation of a wave packet 
in the x-direction compared with the plane wave [exp(lk x)] solutions of 
Chaturvedi (1976). 

We solve Eq. (14) analytically and numerically to determine the 
eigenvalues and eigenfunctions. We follow the techniques developed by 
Ganguli and Bakshi (1982) to solve Eq. (14) as 


M-+ Q(xH - 0, 

3x Z 


[ (m 2 -fl 2 -k 2 c 2 J + iu^jkV/L 2 + i(m/M) v^m-^V^/Ljk 2 


We expand Q(x) around 5 = x-x q * 0 to 0(5 ) where x Q is the value of x for 
which the local growth rate is a maximum, 


( 19 ) 




Now, we have 


—j + (a + n 2 J <Kn) - 0 
3n ' 


where n ■ 5 + Q'/Q" and a ■ Q - Q'/2Q ". The primes indicate the 
o o o o o 

derivatives with respect to x evaluated at x - x Q . 

Equation (20), which is in the form of a Weber equation, has the 
solution (for the l * 0 mode) 


<Kx) ■ <l> exp[- 4 (x - x J 2 /x 2 ] exp [iq(x)J 


where l is the order of the mode and q(x) is a real function. The 
wavepacket peaks at (Ganguli and Bakshi, 1982) 

Mq; /q~i-q~/2j 1/2 j 

- 2 - 2 -- (22 > 

m ° ° Re (-q;'/2j 1/2 

The dispersion relation is given by 


(23) reduces, for the l »■ 0 mode, to 


Equation (23) is anologous to the dispersion relation obtained by Ganguli 
and Bakshi (1982). In Eq. (23) the nonlocal effects are contained in the two 
terms on the right hand side, one of which contains an explicit shear 
dependent term and the other contains a shear independent nonlocal term. 


containing Q' could significantly modify the growth rate. In fact, this 
shear Independent nonlocal term drastically reduces the growth rate in the 
collisionless kinetic case. Whereas, Eq. (24) shows that in the fluid case 
any stabilizing infuelnce by the magnetic shear is due only from the shear 
dependent term whiche indicates that large magnetic shear is probably needed 



to stabilize the instability. In fact, the results given in the next 
section support this observation. Equation (24) is solved to obtain the 
eigenfrequencles and the results are presented in the next section. 


III. RESULTS 

We solve the dispersion equation given by Eq. (24) numerically in 
various parameter domains. These results are verified by directly solving 
Eq. (14) using a numerical shooting code. We present the former results in 
the following. In Fig. 1 we plot the growth rate normalized to the ion 

cyclotron frequency (yAO versus the normalized wavenumber (k p ) for 
i y s 

V^/c g ■ 50, v e /^ e * 0.01, and v^/Q^ * 0.01. In plotting these results we 

treated P g /L g as an external parameter. The solid line represents t' 
shear-free case and the dashed line represents the case where P g /L g “ 0*01 
We see from the figure that inodes with small k^p g are not strongly affec 
by the magnetic shear. The growth rate for k^p g * 0*3 reduces from 0.24 -r 
0.20 for p g /L g “ 0.01, whereas the growth rate for k^p g ■ 0.7 reduces by 
about 30% from 0.75 to 0.55. 

In order to throw some light on the magnetic shear required to 
significantly reduce the growth rate, we plot in Fig. 2, the normalized 
growth rate versus the normalized shear length (p g /L g ), treating p g /L g as 
an external parameter. We choose k p g ■ 0.3, and 0.5 and plot the growth 
rate for the parameters of Fig. 1. This figure shows that moderate to 
strong magnetic shears, p g /L g > 0.01, are required to reduce the growth 
rate significantly. For example, for k p g ■ 0.3 the growth rate is reduced 
by about 50% for p g /L g ■ 0.04 and the mode Is stabilized for p g /L g > 0.065, 

and for k p »0.5 the growth rate is reduced by 50% for p /L = 0.04 but 
y s s s 

the mode is stabilized for P g /L g ■ 0.078. However, the theory begins to 

break down for large shear value P g /L g >0.1. We find that for large 

shears, the growth rate drops linearly as a function of inverse shear 
length. For very small shears, p g /L g < 10 the growth rate remains 
essentially constant and the growth rate decreases noticeably for P g /L g > 
0.01. In an ionospheric environment B = 0.5 Gauss and p 3 2.5 meters, and 

Eqs. (1) and (3) show that p /L ■ 0.01 and 1.0 correspond to currents of 
? v s s 

0.1 A/m and 1 A/m , respectively. 







In Fig. 3 we present the wave eigenfunctions for the following 
-2 -2 

parameters: V, - 50 c , v ./fl - 10 , v./£J. - 10 and for k p - 0.3, 

d seie ii ys 

Y/8^ ■ 0.24. The solid line represents the real part of the wave function 

($) and the dashed line |$|. Two important features are to be noted in 

this figure. One feature is that the wavepacket localizes around x Q /p g « 1, 
i.e., for small shear lengths k^/k^ is « 1, indicating that the mode is 
almost perpendicularly propagating. This can also be seen from Eq. (19) 
which yields 

k 2 /k y ' 2Kk y V d Xp g /L g j « 1 for P g /L g - 10" 7 . 

The second feature is that the width of the wave packet is of the order 
200 p suggesting a localization region of 500 meters for p * 2.5 meters 

3 S 

(corresponding to NO ions in the ionosphere). 

The effect of self-consistent magnetic shear is understood by plotting 

the normalized growth rate versus the normalized drift velocity. For this 

purpose, we use Eq. (4) to express the shear length in terms of the drift 

velocity and solve Eq. (24). In figure 4, we plot the growth rate 

for a - 10 10 \ and 10 * (w^ /u>^ =» 10^, lo\ and 10 4 , respectively; 

pe ce 

c g » 500 m/s) and for v e /^ e = 0.01, * 0.01, and k y P g = 0.3. We also 

give the shear-free local growth rate (curve A) to compare with the growth 
rate with shear. Three points are to be noted in this figure: (1) a » lo""^ 
corresponds to ionospheric parameters. We find that the growth rate is not 
too different from the shear-free case. In fact, the growth rate curve 
overlaps curve A and eventually turns around for V^/cg > 500. (2) Curve B, 
which gives the growth rate for a - 10 shows that the growth rate drops 

significantly beyond drift velocities > 150 c ; and (3) curve C, which 

-4 

represents the growth rate for a *>■ 10 , shows that the optimum drift 

velocity (the drift velocity for which the growth rate is a maximum), 80 c g , 
is much smaller for this case. This leads to the conclusion that the 

optimum drift velocity decreases as the parameter a is increased. At large 

drift velocities, the self-consistent magnetic shear produced is large 

enough that the mode is stabilized and the growth rate tends to zero. 





iv. r. scussion 

We now apply the results to the high latitude Ionospheric E region. We 

choose v ,/n - io“ 2 , v /» « 10 ~ 2 , and V ' 20 c . For V. - 20 c_ the 

11 ee d s a 8 

corresponding shear length Is about 1000 km (N0 + ions, n - 10 5 cm' 3 ). We 

find that the growth rate is 0.9 s' 1 , for modes with wavelengths of 50 

meters (k p * 0.3, and p ■ 2.5 meters), which Is comparable with the local 
y s s 

growth rate. From the plots of the wave functions we see that the mode is 
localized In the north-south direction in a region of Ax ' 200 p g - 500 m. 
Moreover, we find that ky » k x (k^ - 1/Ax) which indicates the strong two 
dimensional structure of the mode in the plane perpendicular to the magnetic 
field during the linear phase of the instability. 

We conclude that the magnetic shear corresponding to normal auroral 

2 

conditions with field aligned currents of the order of pA/ra is small 

(L g - 1000 km) and thus does not apprecibly alter the current driven 

collislonal ion cyclotron mode discussed by Chaturvedi (1976). Even strong 

2 

currents under some disturbed conditions, - 135 pA/m , reported by Burke et 

al. (1983) do not produce significant shear (L g ■* 300 km) in the magnetic 

field to have any effect on the CIC. Currents of the order mA/m 2 or more 

possibly produce significant effects. 

Finally, we discuss the limitations of the present theory. (1) The 

mode structure equation (Eq. 18) and the subsequent dispersion relation (Eq. 

23) were derived under the assumption that k g /ky < 1 based on the premise 

that the ion parallel motion (V»(V^) term In the ion continuity equation] 

may not be important. Thus the results presented do not contain ion 

parallel motion effects. However, we extended the theory to include the ion 

parallel motion and found no significant changes in the way the magnetic 

shear affects the mode. The local growth rate (in the limit p /L - 0 ) 

with the ion parallel motion is smaller by about 10Z; the growth rate with 

magnetic shear is consistently smaller but has similar behaviour as shown 

earlier. (2) For larger shears, p /L '0.1, the parabolic expansion of the 
s s 

potential function (Eq. 23), used to derive the analytical dispersion 
relation, may not be adequate. This Is due to the fact that the wave 
function spreads out and samples non-parabolic part of the potential. 
Numerical solutions of (Eq. 18) do confirm the spreading of the wave 
packet. Furthermore, the wave packet localizes in the region where the 
effective (k g /ky) '1 because magnetic shear introduces an effective k g . (3) 

9 






We have done an MHD analysis and arrived at the mode structure equation with 
a potential term, Q(x) [Eq. (18)], whose derivative vanishes at x Q , whereas, 
In the collisionless case the derivative of Q(x) at x Q is finite, thus 
drastically affecting the collisionless kinetic ion cyclotron instability 
(Ganguli and Bakshi, 1982). To further examine the above three aspects we 
will present the kinetic analysis of the collisonal ion cyclotron waves in 
the presence of magnetic shear in a future paper. For higher drift 
velocities (comparable to the E * B drift velocities) where the self- 
consistent magnetic shear plays a significant role in determining the mode 
structure, a complete treatment is needed which includes an electric field 
along the y-direction and the related velocity shear. 

In conclusion, we have examined the influence of magnetic shear on the 
collisional current driven ion cyclotron instability. Self consistent 
magnetic shear corresponding to moderate drift velocities (V d ^0.5 - 5 
km/s) near the threshold velocity for the instability does not have 
significant effects on the instability. We find that the mode is almost 
perpendicularly propagating and is localized in a region extending upto a 
few hundred kilometers in the direction of the magnetic field under auroral 

conditions. However, in domains where the plaraa density is high such that 
2 2 

the parameter tOp^/o)^ is large, we find that the growth rate of a 
particular mode maximizes at an optimum drift velocity much larger than the 
sound speed. Finally, we find that strong shears (p g /t. s > 0.05) signifi¬ 
cantly reduce the growth rate and stabilize the collisional current driven 
ion cyclotron instability. This strong shear corresponds to large parallel 
drift velocities (V d » c g ), as seen in fig. 4. 
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Figure 1. Plot of normalized growth rate (y/ 0^) versus the normalized 
wave number (k^p s ). The parameters used are v g /Q e * 0.01, 

Vj/0^ - 0.01, and V^/c g « 50. The solid line represents the 
shear free case, and the dashed line represents the case where 








r / Q . i 



Figure 2. Normalized growth rate (yAl^) versus p g /L g treating p g /L g as an 

external parameter for V^/Cg - 50, ly> g ■ 0.3, and 0.5, and for 

v /n - 0.01, and v,/fl. - 0.01. 
e e i i 
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Figure 3. Plot of wave eigenfunctions for the same parameters as Fig. 1 
and for k p ■ 0.3, y/n, - 0.25, and p /L - 0.001. 
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v d /c s 

Figure 4. Growth rate versus the drift velocity curves with and without 
self-consistent magnetic shear for the following parameters: 
k y p s * v e^ n e * 0,01 » and “ 0*01* Curve A represents 

the shear-free case and a ■ 10 Curves B and C are for a ■ 

10" 5 and 10"\ respectively, and these correspond to /w^ ■ 

O A pe ce 

10 3 and 10 . 
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COMMANDER 

U.S. ARMY SATCOM AGENCY 
FT. MONMOUTH, NJ 07703 

OICY ATTN DOCUMENT CONTROL 

COMMANDER 

U.S. ARMY MISSILE INTELLIGENCE AGENCY 
REDSTONE ARSENAL, AL 35809 
OICY ATTN JIM GAMBLE 

DIRECTOR 

U.S. ARMY TRADOC SYSTEMS ANALYSIS ACTIVITY 
WHITE SANDS MISSILE RANGE, NM 88002 
OICY ATTN ATAA-SA 
OICY ATTN TCC/F. PAYAN JR. 

OICY ATTN ATTA-TAC LTC J. HESSE 

COMMANDER 

NAVAL ELECTRONIC SYSTEMS COMMAND 
WASHINGTON, D.C. 20360 

OICY ATTN NAVALEX 034 T. HUGHES 
OICY ATTN PME 117 
OICY ATTN PME 117-T 
OICY ATTN CODE 5011 

COMMANDING OFFICER 
NAVAL INTELLIGENCE SUPPORT CTR 
4301 SUITLAND ROAD, BLDG. 5 
WASHINGTON, D.C. 20390 

OICY ATTN MR. DUBBIN STIC 12 

OICY ATTN NISC-50 

OICY ATTN CODE 5404 J. GALET 

COMMANDER 

NAVAL OCCEAN SYSTEMS CENTER 
SAN DIECO, CA 92152 
OICY ATTN J. FERGUSON 
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SAUL RESEARCH LABORATORY 

WASHINGTON, D.C. 20375 

OICY ATTN CODE 4700 S. L. Ossakov 

26 CYS IF UNCLASS. 1 CY IF CLASS) 
OICY ATTN CODE 4701 I Vltkovitsky 
OICY ATTN CODE 4780 J. Huba (100 
CYS IF UNCLASS, I CY IF CLASS) 
OICY ATTN CODE 7500 

OICY ATTN CODE 7550 

OICY ATTN CODE 7580 

OICY ATTN CODE 7551 

OICY ATTN CODE 7555 

OICY ATTN CODE 4730 E. MCLEAN 

OICY ATTN CODE 4108 

OICY ATTN CODE 4730 B. RIPIN 

20CY ATTN CODE 2628 

COMMANDER 

NAVAL SEA SYSTEMS COMMAND 

WASHINGTON, D.C. 20362 
OICY ATTN CAPT R. PITKIN 

COMMANDER 

NAVAL SPACE SURVEILLANCE SYSTEM 

DAHLGREN, VA 22448 

OICY ATTN CAPT J.H. BURTON 


COMMANDER 

AEROSPACE DEFENSE COMMAND/DC 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 

OICY ATTN DC MR. LONG 

COMMANDER 

AEROSPACE DEFENSE COMMAND/XPD 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 
OICY ATTN XPDQQ 
OICY ATTN XP 


AIR FORCE GEOPHYSICS LABORATORY 
HANSCOM AFB, MA 01731 

OICY ATTN OPR HAROLD GARDNER 
OICY ATTN LKB KENNETH S.W. CHAMPION 
OICY ATTN OPR ALVA T. STAIR 
OICY ATTN PHD JURGEN BUCHAU 
OICY ATTN PHD JOHN P. MULLEN 


AF WEAPONS LABORATORY 
KIRTLAND AFT, NM 87117 
OICY ATTN SUL 

OICY ATTN CA ARTHUR H. GUENTHER 
OICY ATTN NTYCE 1LT. G. KRAJEI 


OFFICER-IN-CHARGE 
NAVAL SURFACE WEAPONS CENTER 
WHITE OAK, SILVER SPRING. MD 20910 
OICY ATTN CODE F31 

DIRECTOR 

STRATEGIC SYSTEMS PROJECT OFFICE 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20376 
OICY ATTN NSP-2141 
OICY ATTN NSSP-2722 FRED WIMBERLY 

COMMANDER 

NAVAL SURFACE WEAPONS CENTER 
DAHLGREN LABORATORY 
DAHLGREN, VA 22448 

OICY ATTN CODE DF-14 R. BUTLER 

OFFICER OF NAVAL RESEARCH 
ARLINGTON, VA 22217 
OICY ATTN CODE 465 
OICY ATTN CODE 461 
OICY ATTN CODE 402 
OICY ATTN CODE 420 
OICY ATTN CODE 421 


AFTAC 

PATRICK AFB, FL 32925 
OICY ATTN TF/MAJ WILEY 
OICY ATTN TN 

AIR FORCE AVIONICS LABORATORY 
WRIGHT-PATTERSON AFB, OH 45433 
OICY ATTN AAD WADE HUNT 
OICY ATTN AAD ALLEN JOHNSON 

DEPUTY CHIEF OF STAFF 
RESEARCH, DEVELOPMENT, & ACQ 
DEPARTMENT OF THE AIR FORCE 
WASHINGTON, D.C. 20330 
OICY ATTN AFRDQ 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 
OICY ATTN J. DEAS 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION/YSEA 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01732 
OICY ATTN YSEA 
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HEADQUARTERS 

ELECTRONIC SYSTEMS D1YISIOS/DC 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 

01CY ATTN DCKC MAJ J.C. CLARK 

COMMANDER 

FOREIGN TECHNOLOGY DIVISION, AFSC 
WRIGHT-PATTERSON AFB, OH 45433 
01CY ATTN NICD LIBRARY 

01CY ATTN ETDP B. BALLARD 

COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
GRIFFISS AFB, NY 13441 

01CY ATTN DOC LIBRARY/TSLD 
01CY ATTN OCSE V. COYNE 

SAMSO/SZ 

POST OFFICE BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 
(SPACE DEFENSE SYSTEMS) 

01CY ATTN SZJ 

STRATEGIC AIR COMMAND/XPFS 
OFFUTT AFB, NB 68113 

01CY ATTN ADWATE MAJ BRUCE BAUER 
01CY ATTN NRT 

01CY ATTN DOK CHIEF SCIENTIST 

SAMSO/SK 
P.O. BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 

01CY ATTN SKA (SPACE COMM SYSTEMS) 
M. CLAVIN 

SAMSO/MN 

NORTON AFB, CA 92409 
(MINUTEMAN) 

01CY ATTN MNNL 

COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
HANSCOM AFB, MA 01731 

01CY ATTN EEP A. LORENTZEN 

DEPARTMENT OF ENERGY 
LIBRARY ROOM G-042 
WASHINGTON, D.C. 20545 

01CY ATTN DOC CON FOR A. LABOWITZ 


DEPARTMENT OF ENERGY 
ALBUQUERQUE OPERATIONS OFFICE 
P.O. BOX 5400 
ALBUQUERQUE, NM 87115 

01CY ATTN DOC CON FOR D. SHERWOOD 

EG&G, INC. 

LOS ALAMOS DIVISION 

P.O. BOX 809 

LOS ALAMOS, NM 85544 

01CY ATTN DOC CON FOR J. BREEDLOVE 

UNIVERSITY OF CALIFORNIA 
LAWRENCE LIVERMORE LABORATORY 
P.O. BOX 808 
LIVERMORE, CA 94550 

01CY ATTN DOC CON FOR TECH INFO DEPT 

01CY ATTN DOC CON FOR L-389 R. OTT 

01CY ATTN DOC CON FOR L-31 R. HAGER 

01CY ATTN DOC CON FOR L-46 F. SEWARD 

LOS ALAMOS NATIONAL LABORATORY 

P.O. BOX 1663 

LOS ALAMOS, NM 87545 

01CY ATTN DOC CON FOR J. WOLCOTT 

01CY ATTN DOC CON FOR R.F. TASCHEK 

01CY ATTN DOC CON FOR E. JONES 

01CY ATTN DOC CON FOR J. MALIK 

01CY ATTN DOC CON FOR R. JEFFRIES 

OICY ATTN DOC CON FOR J. ZINN 

01CY ATTN DOC CON FOR P. KEATON 

OICY ATTN DOC CON FOR D. WESTERVELT 
OICY ATTN D. SAPPENFIELD 

SAND1A LABORATORIES 
P.O. BOX 5800 
ALBUQUERQUE, NM 87115 

OICY ATTN DOC CON FOR W. BROWN 

OICY ATTN DOC CON FOR A. THORNBROUGH 

OICY ATTN DOC CON FOR T. WRIGHT 

OICY ATTN DOC CON FOR D. DAHLGREN 

OICY ATTN DOC CON FOR 3141 

OICY ATTN DOC CON FOR SPACE PROJECT DIV 

SANDIA LABORATORIES 
LIVERMORE LABORATORY 
P.O. BOX 969 
LIVERMORE, CA 94550 

OICY ATTN DOC CON FOR B. MURPHEY 

OICY ATTN DOC CON FOR T. COOK 

OFFICE OF MILITARY APPLICATION 
DEPARTMENT OF ENERGY 
WASHINGTON, D.C. 20545 

OICY ATTN DOC CON DR. YO SONG 



OTHER GOVERNMENT BOEING COMPANY, THE 

P.O. BOX 3707 

DEPARTMENT OF COMMERCE SEATTLE, WA 98124 

NATIONAL BUREAU OF STANDARDS OICY ATTN G. KEISTER 

WASHINGTON, D.C. 20234 01CY ATTN D. MURRAY 

01CY (ALL CORRES: ATTN SEC OFFICER FOR) 01CY ATTN G. HALL 

01CY ATTN J. KENNEY 


INSTITUTE FOR TELECOM SCIENCES 
NATIONAL TELECOMMUNICATIONS & INFO ADMIN 
BOULDER, CO 80303 

01CY ATTN A. JEAN (UNCLASS ONLY) 

01CY ATTN W. UTLAUT 
01CY ATTN D. CROMBIE 
01CY ATTN L. BERRY 

NATIONAL OCEANIC & ATMOSPHERIC ADMIN 
ENVIRONMENTAL RESEARCH LABORATORIES 
DEPARTMENT OF COMMERCE 
BOULDER, CO 80302 
01CY ATTN R. GRUBB 
01CY ATTN AERONOMY LAB G. REID 


DEPARTMENT OF DEFENSE CONTRACTORS 

AEROSPACE CORPORATION 
P.O. BOX 92957 
LOS ANGELES, CA 90009 
01CY ATTN I. GARFUNKEL 
01CY ATTN T. SALMI 
01CY ATTN V. JOSEPHSON 
01CY ATTN S. BOWER 
OICY ATTN 0. OLSEN 

ANALYTICAL SYSTEMS ENGINEERING COR? 
5 OLD CONCORD ROAD 
BURLINGTON, MA 01803 

OICY ATTN RADIO SCIENCES 

AUSTIN RESEARCH ASSOC., INC. 

1901 RUTLAND DRIVE 
AUSTIN, TX 78758 
OICY ATTN L. SLOAN 
OICY ATTN R. THOMPSON 

BERKELEY RESEARCH ASSOCIATES, INC. 
P.O. BOX 983 
BERKELEY, CA 94701 
OICY ATTN J. WORKMAN 
OICY ATTN C. PRETTIE 
OICY ATTN S. BRECHT 


CHARLES STARK DRAPER LABORATORY, INC. 
555 TECHNOLOGY SQUARE 
CAMBRIDGE, MA 02139 
OICY ATTN D.B. COX 
OICY ATTN J.P. GILMORE 

COMSAT LABORATORIES 
LINTHICUM ROAD 
CLARKSBURG, MD 20734 
OICY ATTN G. HYDE 

CORNELL UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 
ITHACA, NY 14850 
OICY ATTN D.T. FARLEY, JR. 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75080 
OICY ATTN H. LOGSTON 
OICY ATTN SECURITY (PAUL PHILLIPS) 

EOS TECHNOLOGIES, INC. 

606 WUshlre Blvd. 

Santa Monica, Calif 90401 
OICY ATTN C.B. GABBARD 

ESL, INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
OICY ATTN J. ROBERTS 
OICY ATTN JAMES MARSHALL 

GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
VALLEY FORGE SPACE CENTER 
GODDARD BLVD KING OF PRUSSIA 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 

OICY ATTN M.H. BORTNER SPACE SCI LAB 

GENERAL ELECTRIC COMPANY 
P.O. BOX 1122 
SYRACUSE, NY 13201 
OICY ATTN F. REIBERT 
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GENERAL ELECTRIC TECH SERVICES CO. INC. 
HMES 

COl’RT STREET 
SYRACUSE, NY 13201 
01CY ATTN G. MILLMAN 

GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 

(ALL CLASS ATTN: SECURITY OFFICER) 
01CY ATTN T.N. DAVIS (UNCLASS ONLY) 
01CY ATTN TECHNICAL LtBRARY 
01CY ATTN NEAL BROWN (UNCLASS ONLY) 

GTE SYLVANIA, INC. 

ELECTRONICS SYSTEMS GRP-EASTERN DIV 
77 A STREET 
NEEDHAM, MA 02194 

01CY ATTN DICK STEINHOF 

HSS, INC. 

2 ALFRED CIRCLE 
BEDFORD, MA 01730 

01CY ATTN DONALD HANSEN 

ILLINOIS, UNIVERSITY OF 
107 COBLE HALL 
150 DAVENPORT HOUSE 
CHAMPAIGN, IL 61820 

(ALL CORRES ATTN DAN MCCLELLAND) 

01CY ATTN K. YEH 

INSTITUTE FOR DEFENSE ANALYSES 
1801 NO. BEAUREGARD STREET 
ALEXANDRIA, VA 22311 
01CY ATTN J.M. AEIN 
01CY ATTN ERNEST BAUER 
01CY ATTN HANS WOLFARD 
01CY ATTN JOEL BENGSTON 

INTL TEL & TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 
01CY ATTN TECHNICAL LIBRARY 

JAYCOR 

11011 TORREYANA ROAD 
P.O. BOX 85154 
SAN DIEGO, CA 92138 
OICY ATTN J.L. SPERLING 


JOHNS HOPKINS UNI VERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAUREL, MD 20810 

OICY ATTN DOCUMENT LIBRARIAN 
OICY ATTN THOMAS POTEMRA 
OICY ATTN JOHN DASSOULAS 

KAMAN SCIENCES CORP 
P.O. BOX 7463 

COLORADO SPRINGS, CO 80933 
OICY ATTN T. MEAGHER 

KAMAN TEMPO -CENTER FOR ADVANCED STUDIES 
816 STATE STREET (P.O DRAWER QQ) 

SANTA BARBARA, CA 93102 
OICY ATTN DASIAC 
OICY ATTN WARREN S. KNAPP 
OICY ATTN WILLIAM MCNAMARA 
OICY ATTN B. GAMBILL 

LINKABIT CORP 
10453 ROSELLE 
SAN DIEGO, CA 92121 
OICY ATTN IRWIN JACOBS 

LOCKHEED MISSILES & SPACE CO., INC 
P.O. BOX 504 
SUNNYVALE, CA 94088 
OICY ATTN DEPT 60-12 
OICY ATTN D.R. CHURCHILL 

LOCKHEED MISSILES & SPACE CO., INC. 

3251 HANOVER STREET 
PALO ALTO, CA 94304 

OICY ATTN MARTIN WALT DEPT 52-12 
OICY ATTN W.L. IMHOF DEPT 52-12 
OICY ATTN RICHARD G. JOHNSON DEPT 52-12 
OICY ATTN J.B. CLADIS DEPT 52-12 

MARTIN MARIETTA CORP 
ORLANDO DIVISION 
P.O. BOX 5837 
ORLANDO, Ft 32805 
OICY ATTN R. HEFFNER 


M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, MA 02173 

OICY ATTN DAVID M. TOWLE 
OICY ATTN L. LOUGHLIN 
OICY ATTN D. CLARK 
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MCDONNEL DOUGLAS CORPORATION 
5301 BOLSA AVENUE 
HUNTINGTON BEACH, CA 92647 
OICY ATTN N. HARRIS 
OICY ATTN J. MOULE 
01CY ATTN GEORGE MROZ 
01CY ATTN W. OLSON 
01CY ATTN R.W. HALPRIN 
01CY ATTN TECHNICAL LIBRARY SERVICES 

MISSION RESEARCH CORPORATION 
735 STATE STREET 
SANTA BARBARA, CA 93101 
01CY ATTN P. FISCHER 
01CY ATTN W.F. CREVIER 
01CY ATTN STEVEN L. GUTSCHE 
01CY ATTN R. BOGUSCH 
01CY ATTN R. HENDRICK 
01CY ATTN RALPH KILB 
01CY ATTN DAVE SOWLE 
01CY ATTN F. FAJEN 
01CY ATTN M. SCHEIBE 
01CY ATTN CONRAD L. LONGMIRE 
01CY ATTN B. WHITE 

MISSION RESEARCH CORP. 

1720 RANDOLPH ROAD, S.E. 

ALBUQUERQUE, NEW MEXICO 87106 
01CY R. STELLINGWERF 
01CY M. ALME 
01CY L. WRIGHT 

MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

01CY ATTN JOHN MORGANSTERN 
01CY ATTN G. HARDING 
01CY ATTN C.E. CALLAHAN 

MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
01CY ATTN W. HALL 
01CY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP 
12340 SANTA MONICA BLVD. 

LOS ANGELES, CA 90025 

01CY ATTN E.C. FIELD, JR. 


PENNSYLVANIA STATE UNIVERSITY 
IONOSPHERE RESEARCH LAB 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 
(NO CLASS TO THIS ADDRESS) 

01CY ATTN IONOSPHERIC RESEARCH LAB 

PHOTOMETRICS, INC. 

4 ARROW DRIVE 
WOBURN, MA 01801 

01CY ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 3027 
BELLEVUE, WA 98009 

01CY ATTN E.J. FREMOUW 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 10367 
OAKLAND, CA 94610 
ATTN A. THOMSON 

RAD ASSOCIATES 
P.O. BOX 9695 
MARINA DEL REY, CA 90291 
01CY ATTN FORREST GILMORE 
01CY ATTN WILLIAM B. WRIGHT, JR. 
01CY ATTN ROBERT F. LELEVIER 
01CY ATTN WILLIAM J. KARZAS 
OICY ATTN H. ORY 
01CY ATTN C. MACDONALD 
OICY ATTN R. TURCO 
OICY ATTN L. DeRAND 
OICY ATTN W. TSAI 

RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
OICY ATTN CULLEN CRAIN 
OICY ATTN ED BEDROZIAN 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 
OICY ATTN BARBARA ADAMS 

RIVERSIDE RESEARCH INSTITUTE 
330 WEST 42nd STREET 
NEW YORK, NY 10036 

OICY ATTN VINCE TRAPANI 
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SCIENCE APPLICATIONS, ISC. 
1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 

01CY ATTN LEWIS M. LINSON 
01CY ATTN DANIEL A. HAMLIN 
01CY ATTN E. FRIEMAN 
01CY ATTN E.A. STRAKER 
01CY ATTN CURTIS A. SMITH 
01CY ATTN JACK MCDOUGALL 

SCIENCE APPLICATIONS, INC 
1710 GOODRIDGE DR. 

MCLEAN, VA 22102 
ATTN: J. COCKAYNE 

SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 94025 

01CY ATTN DONALD NEILSON 
01CY ATTN ALAN BURNS 
01CY ATTN G. SMITH 
01CY ATTN R. TSUNODA 
01CY ATTN DAVID A. JOHNSON 
01CY ATTN WALTER G. CHESNUT 
01CY ATTN CHARLES L. RINO 
01CY ATTN WALTER JAYE 
01CY ATTN J. VICKREY 
01CY ATTN RAY L. LEADABRAND 
01CY ATTN G. CARPENTER 
01CY ATTN G. PRICE 
01CY ATTN R. LIVINGSTON 
01CY ATTN V. GONZALES 
01CY ATTN D. MCDANIEL 

TECHNOLOGY INTERNATIONAL CORP 
75 WIGGINS AVENUE 
BEDFORD, MA 01730 
01CY ATTN W.P. BOQUIST 

TOYON RESEARCH CO. 

P.O. Box 6890 
SANTA BARBARA, CA 93111 
01CY ATTN JOHN ISE, JR. 

01CY ATTN JOEL GARBARINO 


TRW DEFENSE & SPACE SYS GROU 
ONE SPACE PARK 
REDONDO BEACH, CA 90278 
01CY ATTN R. K. PLEBUCH 
01CY ATTN S. ALTSCHULER 
01CY ATTN D. DEE 
01CY ATTN D/ STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH BEDFORD STREET 
BURLINGTON, MASS 01803 
01CY ATTN W. REIDY 
01CY ATTN J. CARPENTER 
01CY ATTN C. HUMPHREY 
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IONOSPHERIC MODELING DISTRIBUTION LIST 
(UNCLASSIFIED ONLY) 

PLEASE DISTRIBUTE ONE COPY TO EACH OF TEE FOLLGVING PEOPLE (UNLESS OTHERWISE • 

NOTED) 


NAVAL RESEARCH LABORATORY 
WASHINGTON, D.C. 20375 

Dr. P. MANGE - CODE 4101 
Dr. P. GOODMAN - CODE 4180 

A.F. GEOPHYSICS LABORATORY 
L.G. HANSCOM FIELD 
BEDFORD, MA 01730 
DR. T. ELKINS 
DR. W. SWIDER 
MRS. R. SAGALYN 
DR. J.M. FORBES 
DR. T.J. KENESHEA 
DR. W. BURKE 
DR. H. CARLSON 
DR. J. JASPERSE 

BOSTON UNIVERSITY 
DEPARTMENT OF ASTRONOMY 
BOSTON, MA 02215 
DR. J. AARONS 

CORNELL UNIVERSITY 
ITHACA, NY 14850 
DR. W.E. SWARTZ 
DR. D. FARLEY 
DR. M. KELLEY 

HARVARD UNIVERSITY 
HARVARD SQUARE 
CAMBRIDGE, MA 02138 
DR. M.B. McELROY 
DR. R. LINDZEN 

INSTITUTE FOR DEFENSE ANALYSIS 
400 ARMY/NAVY DRIVE 
ARLINGTON, VA 22202 
DR. E. BAUER 

MASSACHUSETTS INSTITUTE OF 
TECHNOLOGY 

PLASMA FUSION CENTER 
LIBRARY, NW16-262 
CAMBRIDGE, MA 02139 


NASA 

GODDARD SPACE FLICHT CENTER 
GREEKBELT, MD 20771 
DR. K. MAEDA 
DR. S. CURTIS 
DR. M. DUBIN 

DR. N. MAYNARD - CODE 696 
COMMANDER 

NAVAL AIR SYSTEMS COMMAND 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20360 
DR. T. CZUBA 

COMMANDER 

NAVAL OCEAN SYSTEMS CENTER 
SAN DIEGO, CA 92152 

MR. R. ROSE - CODE 5321 

NOAA 

DIRECTOR OF SPACE AND 
ENVIRONMENTAL LABORATORY 
BOULDER, CO 80302 
DR. A. GLENN JEAN 
DR. G.W. ADAMS 
DR. D.N. ANDERSON 
DR. K. DAVIES 
DR. R.F. DONNELLY 

OFFICE OF NAVAL RESEARCH 
800 NORTH QUINCY STREET 
ARLINGTON, VA 22217 
DR. G. JOINER 

PENNSYLVANIA STATE UNIVERSITY 
UNIVERSITY PARK, PA 16802 
DR. J.S. NISBET 
DR. P.R. ROHRBAUGH 
DR. L.A. CARPENTER 
DR. M. LEE 
DR. R. DIVANY 
DR. P. BENNETT 
DR. F. KLEVANS 

SCIENCE APPLICATIONS, INC. 
1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 
DR. D.A. HAMLIN 
DR. E. FRIEMAN 


P 
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STAFFORD UNIVERSITY 
STANFORD, CA 94305 
DR. ?.M. BANKS 

U.S. ARMY ABERDEEN RESEARCH 
AND DEVELOPMENT CENTER 
BALLISTIC RESEARCH LABORATORY 
ABERDEEN, MD 

DR. J. HEIMERL 

GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 
DR. L.E. LEE 

UNIVERSITY OF CALIFORNIA, 
BERKELEY 

BERKELEY, CA 94720 
DR. M. HUDSON 

UNIVERSITY OF CALIFORNIA 
LOS ALAMOS SCIENTIFIC LABORATORY 
J-10, MS-664 
LOS ALAMOS, NM 87545 
DR. M. PONGRATZ 
DR. D. SIMONS 
DR. G. BARASCH 
DR. L. DUNCAN 
DR. P. BERNHARDT 
DR. S.P. GARY 

UNIVERSITY OF MARYLAND 
COLLECE PARK, MD 20740 
DR. K. PAPADOPOULOS 
DR. E. OTT 

JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAUREL, MD 20810 
DR. R. GREENWALD 
DR. C. MENG 

UNIVERSITY OF PITTSBURGH 
PITTSBURGH, PA 15213 
DR. N. ZABUSKY 
DR. M. BIONDI 
DR. E. OVERMAN 


UNIVERSITY OF TEXAS 
AT DALLAS 

CENTER FOR RESEARCH SCIENCES 
P.O. BOX 688 
RICHARDSON, TX 75080 
DR. R. HEELIS 
DR. W. HANSON 
DR. J.P. McCLURE 

UTAH STATE UNIVERSITY 
4TH AND 8TH STREETS 
LOGAN, UTAH 84322 
DR. R. HARRIS 
DR. K. BAKER 
DR. R. SCHUNK 
DR. J. ST.-MAURICE 

PHYSICAL RESEARCH LABORATORY 
PLASMA PHYSICS PROGRAMME 
AHMEDABAD 380 009 
INDIA 

P.J. PATHAK, LIBRARIAN 

LABORATORY FOR PLASMA AND 
FUSION ENERGY STUDIES 
UNIVERSITY OF MARYLAND 
COLLEGE PARK, MD 20742 
JHAN VARYAN HELLMAN, 
REFERENCE LIBRARIAN 
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